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SUMMARY 

A new capability called the AUTOSPC feature has been incorporated in the 
April 1987 release of NASTRAN. It gives the user the option of automatically 
applying single-point constraints for the purpose of removing potential grid 
(and scalar) point singularities that have not been otherwise already 
constrained out. This paper gives details of the implementation of this 
feature and describes its usage, illustrating it with an example problem. 


INTRODUCTION 

In pre-April 1987 releases of NASTRAN, a table data block called GPST 
(Grid Point Singularity Table) was generated by the EMA (Element Matrix 

Assembler) module in all Rigid Formats. This table contains data on potential 
grid (and scalar) point singularities and was based on an examination of the 
stiffness matrix assembled by the EMA module. Based on the constraints 
specified by the user, as processed by the GP4 (Geometry Processor - Phase 4) 
module, this GPST table was subsequently examined by the GPSP (Grid Point 
Singularity Processor) module in order to identify those singularities that 

have been removed by the application of single-point or multipoint 
constraints. Any singularities not so removed were identified for the user's 
information because of potential problems that might be caused in the 
subsequent decomposition process. 

In pre-April 1987 releases of NASTRAN, there was also a supplementary 
functional module called GPSPC (Constrain Stiffness Matrix Singularities). 
While not a part of any of the Rigid Formats, this module was designed to 

replace the GPSP module and to do not only what GPSP was doing, but, in 

addition, to automatically constrain singularities not otherwise constrained. 
However, this GPSPC module has had serious flaws both in design and 
implementation and has never worked properly as advertised in the NASTRAN 
User's Manual . 

In order to resolve the above problem and have an effective means of 
automatically constraining singularities, the AUTOSPC feature has been 
introduced in the April 1987 release of NASTRAN. Details of its 
implementation and usage are given below, along with an illustrative example 
problem. 
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IMPLEMENTATION OF THE CAPABILITY 


For the sake of efficiency and user convenience, the AUTOSPC feature was 
incorporated into NASTRAN by modifying and expanding the code in the existing 
GP4 module. This is the most logical and elegant way of doing this as this 
module processes all constraint data. There is no reason for a separate 
module to handle the new feature. 

The following points summarize the important details of the 

implementation of the capability. 

1. Deletion of old modules 

The GPSP and GPSPC modules existing heretofore in NASTRAN were 
deleted from the code. 

2. Addition of a new module 

The GPST data block that has until now been generated by the EMA 
module has been based on the stiffness matrix assembled by the EMA 
module and has not included the contributions of general elements. 
This has now been corrected by postponing the computation of the 
GPST data block until after the contributions of the general 
elements have been added to the stiffness matrix assembled by the 
EMA module and by basing the computation of GPST on this total 
stiffness matrix. In order to perform this efficiently, a new 
module called GPSTGEN (GPST GENeration) was developed and 
incorporated into NASTRAN. This new module generates the GPST data 
block using the same method as was used earlier by the EMA module. 

3. Changes to the EMA module 

As a result of the above, the EMA module was considerably simplified 
by the deletion of the code associated with the generation of the 
GPST data block. Also, as a result, the number of output data 

blocks generated by the EMA module was reduced from 2 to 1. 

4. Changes to the GP4 module 

Significant changes were made to the GP4 module. These changes are 
described below. 

a. The GP4 module was modified to process the GPST data block 

generated by the new GPSTGEN module. A new subroutine called 
GP4SP was added to the GP4 module for this purpose. If the 
AUTOSPC feature is not turned on (the default), the GP4 module 
processes the GPST data block in the same way as was done 

earlier by the old GPSP module. However, if the user has 
turned on the AUTOSPC feature, then the GP4 module additionally 
applies single-point constraints to remove potential grid (and 
scalar) point singularities as per the scheme outlined later 
under the usage of the feature. 

b. As a result of the above change, the number of output data 

blocks for the GP4 module was increased from four to five. The 
fifth output data block is suitable for processing by the OFP 
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(Output File Processor) module and contains information about 
potential singularities that have not been removed. 

c. In order to enhance user convenience, a new set, called the 

SAUTO-set, was added to the displacement sets already existing 
and recognized in NASTRAN. This new displacement set consists 
of those degrees of freedom that have been constrained by the 
GP4 module by the automatic application of single-point 

constraints when the AUTOSPC feature has been turned on. This 

new set is a subset of the s-set (single-point constraints). 
Until now, the s-set comprised two subsets, namely, the sg-set 
(permanent constraints) and the sb-set (boundary constraints). 
Under the new scheme, the s-set may also include the SAUTO-set. 
The printouts generated by DIAGs 21 and 22 in the GP4 module 
have been expanded to display the degrees of freedom belonging 
to this new SAUTO-set. (The SAUTO-set in DIAG 21 output is 
referred to as the AUTO SPC displacement set in DIAG 22 

output.) These expanded printouts will be of great help to 

users. 

d. In addition, when automatic single-point constraints are 
applied as above, SPC1 cards are generated and printed out 
(and, on option, also punched out), indicating the degrees of 
freedom that have been automatically constrained. This is 
explained further later under the usage of the feature. 

e. Additional messages were introduced into the GP4 module to 
inform the user of the automatic application of single-point 
constraints when the AUTOSPC feature is turned on. 

5. Resolution of a software problem 

While the new capability was being developed, it was discovered that 
the GPST table was never generated in the HEAT approach, except in 
an extremely rare situation. This was traced to a software problem 
in the EMA module. This problem was resolved when the new GPSTGEN 
module was developed. 

6. Changes to the Rigid Format Data Base 

As a result of the work described above, all of the Rigid Formats in 
the Rigid Format Data Base were revamped as follows: 

* The new GPSTGEN module was incorporated. 

* The old GPSP module was deleted. 

* The EMA and GP4 module DMAP statements were modified as per the 
changes outlined above. 

* Other miscellaneous changes called for by the above modifications 
were made. 

Also, all of the COSMI C-suppl ied DMAP ALTER packages were 
appropriately modified in line with the above changes. 
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USAGE OF THE AUTOSPC FEATURE 


The AUTOSPC feature is exercised by the use of an integer DMAP parameter 
called AUTOSPC, normally specified via a PARAM bulk data card. It is 
available in all Rigid Formats. The default value of 0 for this parameter 
causes the new feature to be not used, thereby not changing the status quo . 
The new feature is exercised by specifying "a non-zero value for thTs 
parameter. Acceptable values are 1, 2, -1 and -2. (Unacceptable values cause 
the singularity processing to be skipped.) The meanings of the various values 
for this parameter are described below. 

AUTOSPC = 0 

The AUTOSPC feature is not used. 

AUTOSPC = 1 

All singularities that are not already: 

(a) removed via single-point constraints, or 

(b) removed via multipoint constraints, or 

(c) specified on SUPORT cards, 

are removed by the automatic application of single-point 
constraints. 

A set of SPC1 cards is generated and printed out for the user's 
information and convenience, indicating the singularities that 
have been automatically removed as above. These SPC1 cards 
have the same SPC set ID as the current subcase (or a SPC set 
ID of 1 if the current subcase has no SPC set). 

AUTOSPC = 2 


There are two possible cases. 

Case 1 There are jto omitted degrees of freedom in the current 
subcase. 

This case is handled in the same way as the AUTOSPC = 1 case. 

SPC1 cards are generated and printed out as in the AUTOSPC = 1 
case. 

Case 2 There are omitted degrees of freedom in the current 
subcase. 

This case is handled in the same way as the AUTOSPC = 1 case, 
but with one important difference as explained below. 

All singularities that are not already: 

(a) removed via single-point constraints, or 

(b) removed via multipoint constraints, or 

(c) specified on SUPORT cards, or 
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are removed by the automatic application of single-point 
constraints, but only if the singularity is part of the o-set 
(omitted set) . 

SPC1 cards are generated and printed out as in the AUTOSPC = 1 
case. 


AUTOSPC = -1 

This case is handled in the same way as the AUTOSPC = 1 case, 
except that the SPC1 cards generated are punched out as well as 
being printed out. 

AUTOSPC = -2 


This case is handled in the same way as the AUTOSPC = 2 case, 
except that the SPC1 cards generated are punched out as well as 
being printed out. 

AUTOSPC < -2 or > 2 


These illegal values cause singularity processing to be skipped 
in the GP4 module. 

As can be seen from the above discussion, when there are no omitted 
degrees of freedom in a subcase, the use of AUTOSPC = 1 1' and AUTOSPC = '2' 
both work the same way, by removing all potential singularities that are not 
already part of the s-set (single-point constraints), m-set (multipoint 
constraints) or r-set (degrees of freedom specified on SUPORT cards). 
However, if there are omitted degrees of freedom in a subcase, then AUTOSPC = 
'2' will cause that singularity to be removed only if it is part of the o-set, 
while AUTOSPC = '1' will cause it to be removed even if it is not part of the 
o-set. This distinction is of relevance and interest in dynamic analyses, 
where the omit feature is frequently used to reduce the size of the analysis 
set. The user thus has a choice in the usage of the AUTOSPC feature in such 
situations. 


COMMENTS ON THE USAGE OF THE AUTOSPC FEATURE 

In pre-April 1987 releases of NASTRAN, the user has been required to 
explicitly remove singularities in many common situations. These include such 
cases as the constraining of grid point rotations in solid elements, the 
constraining of grid point rotations normal to the plate in plate elements and 
the constraining of specific grid point components in axisymmetric elements. 
A GRDSET card has been normally used to explicitly specify such constraints. 
The AUTOSPC feature now provides an alternative means of handling such 
situations in a convenient and elegant manner. This is illustrated in the 
example problem discussed below. 


5 


EXAMPLE PROBLEM 


In order to demonstrate the usage of the AUTOSPC feature, NASTRAN 
Demonstration Problem No. D01-07-1A was selected. This problem performs 
static analysis of a spherical shell using toroidal ring elements (see the 
NASTRAN Demonstration Problem Manual for details). Among other data cards, 
this problem employs a GRDSET bulk data card to constrain grid point 
displacements tangential to the shell surface. 

The standard demo was first run by commenting out the GRDSET card in the 
data deck, thereby introducing singularities into the problem. (In addition, 
in order to reduce the amount of output, the output request was limited to the 
displacement output.) The results of this run are shown in Appendix 1. As 
can be seen, the run terminated due to the singularities introduced into the 
problem. 

The above data was then augmented by adding a PARAM bulk data card with 
AUTOSPC = 1 specified on it. The results of this run are shown in Appendix 2. 
As can be seen, this run ran successfully to completion, giving results that 
match those of the standard demo problem. 

The above runs also illustrate some of the new user messages, as well as 
the expanded output resulting from the use of DIAGs 21 and 22. (Both of the 
runs were made on a post-April 1986 release of NASTRAN, containing the AUTOSPC 
feature discussed in this paper.) 


CONCLUSIONS 

A new capability called the AUTOSPC feature has been incorporated in the 
April 1987 release of NASTRAN. It gives the user the option of automatically 
applying single-point constraints for the purpose of removing potential grid 
(and scalar) point singularities that have not been otherwise already 
constrained out. This paper has presented the details of the implementation 
and usage of this new feature, illustrating it with an example problem. 
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APPENDIX 1 

NASTRAN Demonstration Problem No. 001-07-1A 
Modified to Have Singularities 
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APPENDIX 2 

Modified NASTRAN Demonstration Problem No. D01-07-1A 
of Appendix 1 Run With AUTOSPC = 1 
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